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SECTION  I 


PROGRAM  DESCRIPTION 


1.1  INTRODUCTION 

Litton 's  Guidance  and  Control  Systems  Division  is  involved  in  a 
multi-year  program  in  the  research  and  development  of  a nuclear 
magnetic  resonance  (NMR)  gyro  for  use  in  inertial  navigation 
systems . 

Angular  rate  is  sensed  in  an  NMR  gyro  through  a shift  in  the 
Larmor  frequency  of  an  ensemble  of  precessing  atomic  or  nuclear 
magnetic  moments.  In  the  Litton  approach,  two  different  noble 
gas  nuclear  isotopes  are  used  to  eliminate  explicit  dependence 
on  the  magnetic  field.  The  noble  gas  ensembles  attain  a net 
magnetic  moment  through  spin  exchange  collisions  with  an  optically 
oriented  alkali  metal  vapor.  In  addition,  the  alkali  vapor, 
through  a magnetometer  effect,  detects  the  weak  magnetic  fields 
which  are  generated  by  the  precessing  nuclear  moments  of  the  two 
noble  gases,  and  thus  provides  the  gyro  angular  rate  information. 

The  aim  of  the  physics  research  program  reported  here,  is  to 
develop  a stronger  theoretical  and  empirical  understanding  of 
the  alkali-atom-noble-gas  nuclear  spin  exchange  process,  noble 
gas  spin  relaxation  mechanisms  and  other  related  NMR  phenomena 
in  a noble-gas-alkali-metal  vapor  system. 


li 
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Results  of  the  studies  will  be  relevant  not  only  to  the  develop- 
ment of  an  NMR  gyro  but  also  in  improving  the  basic  understanding 
of  noble  gas  alkali  atomic  interactions  and  will  therefore  be  of 
interest  to  many  investigators  in  related  fields. 


1 . 2 OBJECTIVES  AND  STATUS  OF  THE  CURRENT  PROGRAM 

The  research  tasks  described  in  Sections  II  and  III  are  aimed  at 
providing  a better  understanding  of  spin  dynamic  properties  of  a 
polarized  noble  gas  nuclear  ensemble  in  a noble  gas,  alkali  metal 
vapor  system.  Investigations  into  this  area  of  atomic  physics 
have  relevance  not  only  to  NMR  gyro  research,  but  also  to  other 
possible  areas  of  research  such  as  maser,  magnetometers  and 
polarized  nuclear  targets. 

In  Section  II  a theory  for  the  role  of  the  alkali  atom,  noble 
gas  atom,  Van  der  Waals  molecules  in  the  alkali  electronic  noble 
gas  nuclear  spin  exchange  process  is  developed.  Preliminary 
data  is  presented  and  results  are  interpreted  which  display  a 
significant  contribution  to  the  spin  exchange  interaction  from 
the  formation  of  the  alkali-atom-noble-gas-atom  diatomic 
molecules . 


In  Section  III,  the  relaxation  of  polarized  noble  gas  nuclei 
through  the  interaction  of  the  nuclear  quadrupole  moment  with 
electric  field  gradients  is  considered.  A model  for  this  inter- 
action is  presented  which  qualitatively  explains  various  phenomena 

83 

associated  with  the  spin  relaxation  of  Kr 
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The  work  described  in  both  Section  II  and  III  is  part  of  a 
continuing  effort.  A more  quantitative  understanding  of  the 
associated  phenomena  is  the  aim  of  future  work  in  these  areas. 


1.3  PROFESSIONAL  PERSONNEL  ASSOCIATED  WITH  THE  PROGRAM 
The  principal  investigator  of  these  research  efforts  has  been 
Dr.  Charles  H.  Volk.  This  work  has  also  been  strongly  supported 
by  Dr.  John  G.  Mark,  Mr.  Bruce  C.  Grover  and  Dr.  Tae  M.  Kwon. 

We  also  wish  to  acknowledge  the  contributions  of  Mr.  Howard  E. 
Williams  for  the  design,  and  Mr.  Roger  L.  Meyer  for  the  con- 
struction of  various  components  of  the  apparatus  used  in  this 
research  effort. 
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SECTION  II 

EFFECTS  OF  ALKALI-ATOM-NOBLE-GAS-ATOM  MOLECULES 
ON  THE  SPIN  EXCHANGE  PROCESS 
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2.1  INTRODUCTION 

The  spin  exchange  process  between  an  alkali  valence  electron  and 
a noble  gas  nucleus  has  been  explained  in  terms  of  the  contact- 
hyperfine  interaction  between  an  alkali  atom  and  a noble  gas 
nucleus.  The  effective  hamiltonian  governing  this  interaction 
has  been  shown  to  be : ^ 


Hi  = yi  • S 


where  y is  the  strength  of  the  spin  exchange  interaction,  I is 
the  spin  angular  momentum  of  the  noble  gas  nucleus  and  S is  the 
spin  angular  momentum  of  the  alkali  valence  electron.  The  equa- 
tions of  motion  for  observables  of  interest  have  been  derived 

under  H^  for  alkali  and  noble  gas  atoms  interacting  in  binary 
2 3 

type  collisions,  ’ and  the  experimental  measurements  were  found 

3 4 

to  be  in  very  good  qualitative  agreement.  ’ 

Recent  measurements,  however,  of  the  alkali  atom,  noble  gas 
nucleus,  spin  exchange  cross  section  as  a function  of  buffer 
gas  pressure  have  demonstrated  very  strong  effects  which  cannot 
be  explained  by  a theory  which  contains  binary  collision  effects 
only.  These  measurements  indicate  a significant  role  for  the 
alkali-atom-noble-gas-atom  Van  der  Waals  molecule  in  the  spin 
exchange  process. 
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We  will  develop  in  this  report  the  formalism  to  analyze  the  role 

of  the  Van  der  Waals  molecules  in  our  spin  exchange  measurements. 

We  will  adopt  some  simplifying  assumptions,  as  was  done  in  a 
3 

previous  report  , in  order  to  investigate  the  basics  of  the 

molecular  interaction.  Calculations  are  made  for  a system  con- 

2 

sisting  of  an  atom  with  a single  electron  in  a si/2  state  an<3 
no  nuclear  spin  (the  alkali  atom)  and  an  atom  with  no  electronic 
spin,  S = 0,  and  nuclear  spin  1 = 1/2  (the  noble  gas  atom). 


2.2  THE  ALKALI -ATOM- NOBLE-GAS- ATOM  DIATOMIC  MOLECULE 


The  alkali-atom-noble-gas-atom  Van  der  Waals  molecule  has  been 

described  by  several  authors  ° and  recently  its  contributions 

to  the  relaxation  of  spin  polarized  alkali  atoms  have  been 
9 10 

discussed.  ’ For  the  purposes  of  this  report,  we  shall  review 
the  salient  features  of  this  diatomic  system. 


The  alkali-atom- noble-gas-atom  two-body  problem  reduces  to  a 
consideration  of  one  body,  with  reduced  mass: 

n = maV(ma  + mn)  (2) 


where  M^  is  the  mass  of  the  alkali  atom  and  M^^  is  the  mass  of 

2 2 

the  noble  gas  atom.  The  pair  has  energy,  E = h k /2[i,  where  k 
is  the  momentum  of  the  reduced  particle  in  the  effective  potential: 


UN(r,Q)  = U(r,Q)  + N ( N + l)h2/2(i  r2 


2-2 


(3) 


L 
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where  U(r,Q)  is  the  interatomic  potential,  r is  the  internuclear 

separation,  N is  the  relative  orbital  angular  momentum  quant’im 

number  of  the  pair,  and  Q is  the  solid  angle.  Low  energy  scatter- 

ing  experiments  ’ have  shown,  for  alkali  atom,  noble  gas  atom 

pairs  investigated,  that  the  interatomic  potential  is  a central 

6 12 

potential  function:  moreover,  the  Lennard- Jones0  potential 
with  two  adjustable  parameters,  well  depth  and  equilibrium 
separation,  fits  the  experimental  scattering  data  very  well. 

The  potential  function  is  of  the  form: 

U(r)  = €o  [(rn/r)12  " 2(rn/r)6J  (4) 

where  r is  the  value  of  r at  which  U(r)  is  equal  to  its  minimum 
m ^ 

value,  f . 

o 

A typical  plot  of  energy  versus  internuclear  separation  is  shown 
in  Figure  2-1,  for  the  effective  potential  UN(r).  With  reference 
to  Figure  2-1,  it  is  seen  that  the  potential  UN(r)  can  be  divided 
into  distinct  areas  labeled  I,  II,  and  III.  Area  I is  defined 
such  that  UN  < 0.  The  systems  in  this  region  correspond  to  bound 
molecules,  which  are  formed  only  in  three  body  collisions  and 
are  destroyed  in  subsequent  collisions  of  the  complex  with  any 
other  atom.  Region  II  corresponds  to  quasibound  molecules. 
Molecules  in  this  region  have  energy  levels  such  that  0 s E < E(N), 
and  they  are  trapped  by  the  centrifugal  barrier.  Since  the 
energy  levels  are  positive,  quasibound  molecules  can  be  formed 
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in  both  two  and  three  body  collisions.  In  addition  to  the 
collisional  destruction  mode,  quasibound  molecules  are  char- 
acterized by  a finite  lifetime  via  quantum  mechanical  tunneling 
through  the  centrifugal  barrier.  If  the  half-life  for  decay 
through  the  barrier  is  greater  than  the  average  time  between 
collisions  then  the  quasibound  molecules  behave  like  bound 
states.  When  the  time  between  collisions  is  long  compared  to 
the  half-life  for  tunneling,  then  one  must  consider  the  phenomena 
of  long-lived  binary  collisions  whose  properties  differ  from 
those  of  the  sudden  binary  collisions.  Finally,  states  in 
Region  III  are  those  of  unbound  alkali  atom,  noble  gas  atom 
pairs . 

Theoretical  determinations  of  potential  parameters,  such  as  the 
well  depth  and  the  equilibrium  separation,  are  of  two  types: 
dispersion  force  calculations  and  pseudopotential  calculations. 
The  dispersion  potential  yields  only  information  on  the  attrac- 
tive nature  of  the  internuclear  potential  and  is  of  the  form: 

U(r)  = - C6/r6  - Cg/r8  (5) 

Cg  and  Cg  have  been  calculated  for  several  alkali  atom,  noble 
13 

gas  atom  pairs.  The  dispersion  force  calculation  is  of  limited 
value  since  no  estimate  on  the  minimum  value  of  r for  which  the 
potential  remains  valid  can  be  made. 
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Pseudopotential  calculations  have  been  done  by  Pascale  and 
Vanderplanque^  and  Baylis.^  The  Baylis  potential  is  such  that 
it  agrees  with  the  dispersion  force  calculations  at  long  range 
and  with  experimental  determinations,  in  those  cases  for  which 
these  exist,  near  the  bottom  of  the  well.  Table  2-1  lists  some 
relevant  molecular  parameters  based  on  the  Baylis  potential. 

2.3  EFFECTS  OF  THE  MOLECULAR  LIFETIME  ON  THE  EVOLUTION 
OF  THE  SYSTEM 


We  consider  a system  of  alkali  and  noble  gas  atoms  together  with 
one  or  more  buffer  gases,  typically  Nj  and/or  He,  in  the  presence 
of  some  external  magnetic  field,  which  is  used  to  define  a quanti- 
zation axis  for  this  system. ^ The  interaction  between  the  alkali 
and  noble  gas  atoms  is  assumed  to  be  stochastic  and  governed  by 
the  hamiltonian  defined  in  Eq.  (1).  Consider  then  the  evolution 
of  the  wavef unction  of  the  alkali  atom,  noble  gas  atom  system 


from  time  t to  time  t + At.  We  divide  the  diatomic  pairs  into 

two  classes:  ’ (a)  a fraction,  At/T^.  that  interact  at  an 

instant,  t , where  t is  such  that  t < t < t + At  and  T”"'-  is 
o o o f 

the  rate  of  formation  of  the  alkali-atom-noble-gas-atom  inter- 
acting pairs,  and  (b)  a fraction  (1  - At/T^)  which  do  not  evolve 
during  the  same  time  interval.  We  impose  the  following  restric- 
tions At  « Tf  and  At  » T,  where  T is  the  mean  duration  time  of 
the  interaction.  The  first  relationship  demands  that  the  alkali 
atom,  noble  gas  atom  pair  formation  occurs  in  a time  short  com- 
pared to  the  time  over  which  we  look  for  variations  in  the 
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TABLE  2-1.  SOME  RELEVANT  MOLECULAR  PARAMETERS  BASED  ON 

THE  BAYLIS  POTENTIAL3 


Rb-Kr 

Rb-Xe 

1.46 

2.02 

5.08 

5.25 

59,604 


92,346 


142,832 

92 

221,768 


r/nA  b x 10 

/Vt  x 10 


305°K 


The  parameters  we  have  listed  in  Table  2-1  are  described  below: 

eQ  and  r are  the  well  depth  and  equilibrium  separation  for 
the  alkali-noble-gas  potential.  ^ 


nOB  is  the  number  of  rotational  states  of  the  bound  molecule. 

Nmaxg  is  the  maximum  angular  momentum  for  the  bound  molecule; 
nOT  and  NmaxT  are  similarly  defined  for  the  total  states  of 
the  molecule,  bound  plus  quasibound. 

(nAN/nzJ  is  the  concentration  of  molecules  per  alkali  con- 
centration per  Torr  of  noble  gas.  The  subscripts  again 
denote  the  bound  and  total  molecule.  (The  numbers  in  this 
tabJ e were  calculated  with  the  aid  of  References  7 and  8.) 
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wavef unction,  while  the  second  relation  assumes  that  the 
wavefunction  changes  smoothly  in  time,  since  our  observation 
is  made  over  a time  long  compared  to  the  duration  of  any 
interaction. 

The  time  sequence  may  be  pictured  in  the  following  way: 


t 


t t + T t + At 

o o 


In  the  time  intervals,  t - t and  t + At  - t - T,  the  system  is 
governed  by  the  'free'  hamiltonian,  HQ > namely: 


H 


o 


w S 
s z 


+ WTI 

I z 


(6) 


where  wg  and  w^.  are  the  respective  Larmor  frequencies  and  h has 
been  set  to  unity.  In  the  interval  T the  systems  evolve  under 
the  full  hamiltonian: 


H = H + H, 

o 1 

where  H-^  is  defined  in  Eq.  (1). 


(7) 


The  wavefunction  at  time  t + At  is  then  given  in  terms  of  the 
wavefunction  at  time  t: 

-iH  (t+At-tQ-T)  -iHT  -iH  (t  -t) 

e e 


*(t+At) 


e 


*(t) 


(8) 
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It  is  convenient  at  this  point  to  transform  the  wavefunctions 
to  the  'interaction  picture',  Where: 


'Ht+At)  = e 


iH  (t+At) 


*(t+At) 


iH  t 

e ° *(t)  = *(t) 


iH0(t+At) 


Then  by  multiplying  Eq.  (8)  through  on  both  sides  by  e 
and  using  Eq.  (10),  we  find: 


iH  t iH  T -iHT  -iH  t 
*(t+At)  = e ° ° e ° e e °°  *(t) 


Eq.  (11)  then  defines  the  evolution  operator  in  the  1 interaction 
picture 1 : 


iH  t iHT  -iHT  -iH  t 

rJ  v O O O O O 

U(tQ,T,H1)  = e e e e 


If  we  designate  the  density  matrix  for  the  alkali  atom,  noble 
gas  nucleus  system  as  R,  then  in  the  'interaction  picture'  the 
density  matrix  at  time  t+At  is  related  to  the  density  matrix  at 

4-  v.  2.3.9 

time  t by: 


R(t+At)  = ( 1-At/Tf ) R ( t ) 


+ At/T.<  U(t  ,T,H. ) R ( t ) U (t  ,T  H.)>. 

L O J-  O j J.  c 


! 

* 


LLJ  GUIDANCE  & CONTROL  SYSTEMS 

Litton  5500  Canoga  Avenue.  Woodland  Hills.  California  91365 


FSCM  06481 


where  the  first  term  on  the  right-hand  side  of  Eq.  (13)  represents 
the  probability  that  the  system  remains  unchanged  during  the  time 
interval  At.  The  subscript  ' ave ' in  the  brackets  indicates  an 
average  is  to  be  performed  over  the  stochastic  variables  tQ , T and 
H-^.  We  assume  that  all  values  of  tQ  in  the  interval  At  are  equi- 
probable  and  we  thus  choose  a particular  value,  that  is  we  are 
able  to  rotate  the  evolution  operator  to  form  V(T,H^): 


V(T,Hl)  = e 


iH  T -iHT 


Making  the  approximation: 


dR(  t ) 
dt 


R( t+At ) - R(T) 
At 


We  can  now  write  the  differential  equation  for  the  density 
matrix  of  the  alkali  atom,  noble  gas  nucleus  system  in  the 
interaction: 


= - ( 1/Tf ) ^R(  t ) - <V(T,H1)  R ( t ) vt(T.H1)>ave]  U< 

2.3.1  The  Spin  Exchange  Equations 

2 3 

The  density  matrix  of  the  system  is,  as  previously  defined,  ’ 
the  product  of  the  density  matrices,  P A (x)  ^NG ’ two 

subsystems,  i.e.,  those  of  the  alkali  atom  and  the  noble  gas 
nucleus,  respectively.  We  will  designate  the  basis  states 
jmI,ms>  numerically  as: 


2-10 


and  thus  the  elements  of  the  density  matrix,  R,  can  be  expressed 
without  cumbersome  notation. 


tr 


*1 

I 

I 


We  have  assumed  that  the  density  matrix  of  the  system  can  be 
represented  in  terms  of  the  products  of  the  basis  states  of  the 
subsystems,  since  the  coupling  between  the  alkali  and  noble  gas 
atoms  in  the  Van  der  Waals  molecule  is  weak  (see  Table  2-1). 

This  is  not  a bad  approximation  as  long  as  the  wavefunction  over- 
lap effects  are  taken  into  account.1 

We  define  the  following: 


and 


A = e‘iY  T/4 


_ iDT  / 2V  -iw+T  , 2 -iw  T 

B = e Icos  X e + sm  X e 


C = elDT(cos  X)(si 


n X)  ( 


-iw+T  -iw  T 
e - e 


• +„ 


_ -iDT,  , . v.  / -iw'T  -iw  T 

E = e (cos  X) (sin  X)  le  - e 


) 
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(17a) 

(17b) 

(17c) 

( 17d) 


-iDT  / 2V  -iw  T 2 -iw+T\  . 

= e Icos  X e + sin  X e J (17e) 
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where 


D = (w];  - wg  )/2 


(18a) 


w1  = y/4  ± Vy2/4  + °2/4 


(18b) 


tan  X = D +Vr  /4  + 


D2/  4 /{ 


(y/2) 


(18c) 


The  time  evolution  part  of  the  density  matrix,  VRV  can  then  be 
written  as  shown  in  Figure  2-2. 


The  observables  of  interest  for  the  noble  gas  nucleus  are  <1  ^ 
and  ^I  + ^,  which  are  the  expectation  values  of  the  longitudinal 
and  transverse  polarization  of  the  noble  gas  nucleus,  respec- 
tively. These  are  given  in  terms  of  the  density  matrix: 


<1  > = Tr 


M = T ( 


~ 2 lP  11  P 11 
NG  NG  - 2 - 2 


WM 


<I  + > = Tr  ("i+r1  = p 

L J NR  •=•-■=■ 


2 2 


where  I and  I are  the  corresponding  operators. 
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In  the  limit  of  low  magnetic  field,  which  is  a reasonable  limit 
for  the  noble  gas  nucleus  since  the  Larmor  frequencies  are  much 
lower  than  the  spin  exchange  energy  for  the  magnetic  fields  at 
which  we  work,  we  find  the  spin  exchange  equations  from  Eqs.  (16), 
(19)  and  (20): 


<i2>  * -<l/Tf)[<Iz>  - <SZ>]  sin2(TT/2)ave  (21) 

and 

<I  + > = -(l/Tf)[<I  + > - <S+>]  sin2(yT/2)ave 

-i[<Sz><I+>  - <Iz>  <S+>]  sin (y T ) ave  (22) 


We  average  over  the  collision  duration  by  assuming  that  the 
primary  effect  of  the  molecular  formation  on  the  spin  exchange 
process  is  that  of  extending  the  life  of  the  interaction  only. 
Eqs.  (21)  and  (22)  become  after  averaging: 


2 2 

<iz>  = -d/Tf>  7X?  [<Iz>  ' <S^] 

<i+> = - ( i/Tf ) -y—T-2—  r<i+>  - <s+>i 

i + L J 


(23) 


+ (i/Tf>  t JTt2t2  [<Sz>  <I  + > - <h>  <S+>]  (24) 

where  r represents  the  mean  life  of  the  interaction  state.  For 
the  case  of  low  external  magnetic  field,  Eqs.  (23)  and  (24)  have 
exactly  the  same  form  as  for  the  binary  collision  case. 
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2.3.2  The  Exchange  Rate 

From  Eqs.  (23)  and  (24)  we  define  the  exchange  rate  for  the 
exchange  of  spin  angular  momentum  from  the  alkali  atom  to  the 
noble  gas  nucleus: 


-1 

ex 


y2r2 


2(1  + y 2 t2 ) 


(25) 


In  addition  a spin  exchange  cross  section  can  be  defined  as: 


-1 

ex 


N. 


J V . 
ex  rel 


(26) 


where  Na  is  the  alkali  number  density,  o is  the  exchange  cross 
section  and  Vrel  is  the  relative  velocity  between  the  alkali  and 
noble  gas  atoms. 

Eg.  (25)  is  true,  however;  only  if  a single  exchange  process  is 
present,  e.g.,  binary  collisions.  In  general  one  would  expect 
three  types  of  processes,  so  that  the  total  exchange  rate  could 
be  written  as: 


T 1 = T-1  + T 1 + T-1  (27) 

eX  exbound  exquasi-bound  exbinary 

where  the  first  term  in  Eg.  (27)  represents  the  exchange  rate 
through  bound  molecular  complexes,  the  second  term  is  the  con- 
tribution of  the  guasi-bound  molecules  to  the  exchange  process 
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and  the  third  term  is  spin  exchange  due  to  binary,  alkali  atom, 
noble  gas  atom  collisions. 

We  now  consider  each  term  as  a function  of  the  buffer  gas  pres- 
sure, with  the  form  of  the  exchange  rate  for  each  process  being 
given  by  Eq.  (25): 

a.  Bound  state  molecules  ( alkali-atom-noble-gas-atom 

which  exist  in  Region  I of  Figure  2-1):  As  was  dis- 
cussed in  Section  2.2,  the  bound  complexes  are  formed 
in  three-body  collisions  only,  as  a result  of  conserva- 
tion of  energy  and  thus  for  a given  alkali  and  noble 
gas  atom  concentration,  the  formation  rate,  T"1,  is 
proportional  to  the  buffer  gas  density.  Since  the 
binding  energy  of  these  Van  der  Waals  molecules  is  less 
than  a kT,  the  molecular  state  exists  only  until  its 
next  collision  with  any  atom,  which,  due  to  the  con- 
centration of  constituents  in  our  cells,  makes  it 
inversely  proportional  to  the  buffer  gas  density.  The 
overall  effect  then,  for  sufficiently  low  noble  gas 
densities,  is  that  at  low  buffer  gas  densities  the 
exchange  rate  is  proportional  to  the  buffer  gas  density, 
and  for  high  buffer  gas  densities,  is  inversely  propor- 
tional to  the  buffer  gas  density. 
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Quasi-bound  molecules  (complexes  in  Region  II  of 
Figure  2-1):  The  quasi-bound  molecules  are  formed  in 
two  ways:  (1)  tr rough  two-body  alkali  noble  gas  atom 
resonant  collisions  in  which  tunneling  through  the 
centrifugal  barrier  takes  place,  and  (2)  in  three-body 
collisions,  in  the  same  manner  as  for  the  formation  of 
the  bound  state  molecules.  In  case  '1'  the  formation 
rate  is  constant  in  buffer  gas  density  while  in  case  '2' 
the  formation  rate  is  proportional  to  the  buffer  gas 


f 


density.  The  destruction  of  the  quasi-bound  molecules 
can  also  occur  in  two  ways:  (1)  decay  of  the  resonant 
state  via  tunneling  back  through  the  barrier,  and  (2)  in 
a collision  of  the  quasi-bound  complex  with  another 
atom  as  in  the  case  of  the  bound  molecules.  The  dif- 
ferent creation  and  destruction  modes  for  the  quasi- 
bound molecules  result  in  four  different  types  of  buffer 
gas  density  dependencies  in  the  exchange  rate.  However, 
two  of  these  dependencies  are  indistinguishable  from 
the  other  two  exchange  processes,  that  is,  a quasi- 
bound molecule  formed  in  a three-body  collision  and 
then  destroyed  in  a subsequent  collision  with  another 
atom  behaves  like  a bound  molecule,  and  a quasi-bound 
molecule  formed  through  tunneling  which  then  decays 
through  tunneling  has  the  same  density  dependence  as  a 
binary  collision.  In  addition,  previous  experimental 


-» 
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work  has  suggested  that  the  quasi-bound  molecules 
involving  an  alkali  and  one  of  the  heavier  noble  gases 
(Ar,  Kr,  Xe)  do  not  contribute  significantly  in  alkali, 
noble  gas  interactions.  We  assume  then  that  any  buffer 
gas  density  dependence  of  the  exchange  rate  due  to  the 
quasi-bound  molecules  will  be  very  weak, 
c.  Alkali  atom,  noble  gas  atom,  binary  collisions.  This 
process  leads  to  an  exchange  rate  that  is  independent 
of  buffer  gas  density. 

Considering  then  the  exchange  rate  through  only  bound  state 
molecules  and  binary  collisions,  we  parameterize  the  exchange 
rate  as  a function  of  buffer  gas  density,  in  the  limit  of  low 
noble  gas  densities,  as: 


abnb 

n?  + b 
b 


c 


(28) 


where  a,  b and  c are  constants  and  n^  is  the  density  of  the 
buffer  gas.  Examining  Eq.  (28)  then,  we  would  expect  the 
exchange  rate  to  increase  linearly  in  n^  for  low  buffer  gas 
densities  and  then  decrease  as  l/n.^  to  the  constant  binary  rate 
in  the  limit  of  high  buffer  gas  densities,  if  the  alkali  atom, 
noble  gas  nucleus,  spin  exchange  rate  is  influenced  by  the 
presence  of  the  Van  der  Waals  molecules. 
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2 . 4 EXPERIMENTAL  RESULTS 

The  experimental  apparatus  and  data  acquisition  procedure  have 

3 

been  described  previously.  In  brief,  the  Rb  cell  is  mechanized 

1 8 

as  a magnetometer  as  described  by  Cohen-Tannoud ji , et  al, 
and  thus  the  precession  of  the  noble  gas  nuclei  is  detected  as 
an  oscillating  magnetic  field. 

Sample  data  taken  in  this  manner  is  displayed  in  Figure  2-3. 

This  particular  data  was  from  a cell  containing  natural  Rb, 

15  Terr  Nj  > 30  Torr  He  and  0.5  Torr  Xe^^^.  The  horizontal  axis 
on  all  three  plots  represents  time.  The  full  length  of  each  run 
was  approximately  500  seconds.  The  vertical  axis  on  each  plot 
is  the  amplitude  of  the  signal.  Note  that  the  middle  plot  has 
a scale  larger  than  the  other  two  by  a factor  of  two.  The  top 
plot  in  Figure  2-3  represents  data  taken  at  38°C,  the  middle 
plot  that  taken  at  58°C  and  the  bottom  plot  that  taken  at  74.5°C. 

2.4.1  Data  Reduction 

We  typically  perform  two  types  of  experiments:  one  in  which  we 
examine  the  longitudinal  (T^)  pumping  of  the  noble  gas  nuclei 
and  the  other  in  which  we  observe  the  transverse  (T2)  decay  of 
the  nuclear  spin  system.  Longitudinal  effects,  which  are  rep- 
resented by  Eq.  (23),  are  measured  by  pumping  the  nuclear  spin 
system,  through  spin  exchange,  for  various  specified  lengths  of 
time  and  then  recording  the  amplitude  of  the  nuclear  magnetization. 
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Assuming  that  the  nuclear  magnetization  obeys  a simple  exponential 
pumping  law: 


M(t)  = M (1  - 
o 


(29) 


one  can  fit  the  experimental  data  at  a particular  temperature 

for  the  'pump  rate1,  T-"*-;  represented  in  Figure  2-4.  Repeating 

P 

this  as  a function  of  cell  temperature,  and  assuming  the  dominant 
temperature  dependence  of  the  pump  rate  is  due  to  the  spin 
exchange  process,  with  the  temperature  dependence  given  by 
Eg.  (26),  we  have: 


T 1 ( T)  = T"1  + T?1  (T) 
p 1 lex' 


(30) 


and  thus  by  performing  a linear  fit  to  the  pump  rate  data  as  a 
function  of  N^/Vre2  > one  can  find  a value  for  T“i , the  longi- 
tudinal relaxation  rate,  and  the  spin  exchange  cross  section  for 
any  temperature,  as  shown  in  Figure  2-5.  We  note  here  that  we 

fit  the  pump  rate  data  as  a function  of  N./V  . and  not  N„V  . , 

A rel  A rel 

as  might  be  expected  from  Eg.  (26),  since  by  Ref.  1 the  cross 


section  for  spin  exchange  is  proportional  to  V 


-2 
rel  * 


In  our  particular  experiments,  the  alkali  is  Rb,  and  we  use  the 

19 

saturated  rubidium  density  formula  for  the  number  density: 


log  NRb  = -4560/T  - 2.45  log(T)  + 30.98 


(31) 
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TIME  (MIN) 


INITIAL  AMPLITUDE  DATA  IS  FIT  TO  THE  FUNTIONAL  FORM  GIVEN  IN 
EQ.  (29),  THIS  PLOT  REPRESENTS  DATA  TAKEN  FROM  A CELL  CONTAIN- 
ING NATURAL  Rb,  15  TORR  N2.  lOOTorr  He  AND  0.5  Torr  Xe(129)  AT  A 

TEMPERATURE  OF  49.2°C.  A LEAST  SQUARES  FIT  OF  THE  DATA  YIELDED 
A PUMP  TIME,  Tp  = 1727,  THE  ERROR  BARS  REPRESENT  A 5 PERCENT 
i ' UNCERTAINTY  IN  THE  AMPLITUDE. 


. t 


Figure  2-4.  Initial  Noble  Gas  Signal  Amplitude  Versus 

Pump  Time 


Li 
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NRbVRel<10  cm  sec  1 

INTERCEPT  REPRESENTS  THE  TEMPERATURE  INDEPENDENT  DECAY  RATE  AND 
THE  SLOPE  YIELDS  INFORMATION  ON  THE  SPIN  EXCHANGE  CROSS  SECTION. 
ERROR  BARS  REPRESENT  A 5 PERCENT  UNCERTAINTY  IN  THE  DECAY  RATE. 
DATA  IS  TAKEN  FROM  A CELL  CONTAINING  Rb<87),  20.5  Torr  N2  AND  0.05  Torr 
Xe(129K  LEAST  SQUARES  FIT  YIELDED  A TEMPERATURE  INDEPENDENT  RATE 
OF  4.5  x 10-<*  sec-^  AND  oexV^re|  = 6.1  x 10~^  cm4  ^,.-2 


AND  oc 


Figure  2-5.  Noble  Gas  Decay  Rate  Versus  Rb  Density 
Times  Relative  Velocity 
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where  T is  the  cell  temperature  in  degrees  Kelvin.  The 
transverse  decay  rate  of  the  nuclear  spin  system  is  found  by 
fitting  the  data  such  as  that  displayed  in  Figure  2-3  to  a 
decaying  exponential  function.  Referring  to  Eq.  (24),  one  finds 
that  if  the  T ^ of  the  alkali  spin  system  is  much  shorter  than 
that  of  the  nuclear  spin  system,  then  the  decay  of  the  nuclear 
levels  and  the  alkali  electronic  levels  effectively  decouple 
and  Eq.  (24)  becomes: 

<i  + > = -d/Tex*  <I  + > + ^ i/Tex)  <szXI+>  (32) 


The  real  part  of  Eq.  (32)  represents  the  transverse  decay  of  the 
nuclear  levels.  Assuming  that  other  contributions  to  the  trans- 
verse rate  can  be  represented  by  a single  decay  rate,  we  find 
that  the  nuclear  spin  decay  to  be  given  by: 


<ix>-  - [Tex  + (T2rl]  <k> 


where  T2'*'  = T ^ + (TJ,)  ^ is  the  total  transverse  decay  rate  and 
(T^)-'*'  is  the  temperature  independent  component  of  the  total 
transverse  rate.  We  again  assume,  as  was  done  for  the  longitudi- 
nal pumping,  that  the  dominant  temperature  dependence  in  the 
decay  rate  is  due  to  the  exchange  term.  Performing  a linear  fit 
of  T"1  as  a function  of  NA/Vre^  yields  the  temperature  indepen- 
dent T2  rate,  T^,  and  the  transverse  exchange  cross  section. 
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By  Eqs.  (23)  and  (24)  the  exchange  cross  section,  found  by 
analyzing  the  longitudinal  pumping,  is  equal  to  that  found  by 
the  transverse  decay  rate. 

2.4.2  Analysis  of  the  Molecular  Effects 

We  have  initiated  spin  relaxation  measurements  of  Xe-*-^3  in  a 
series  of  cells  designed  to  yield  information  concerning  the 
role  of  the  Rb-Xe  diatomic  molecules  in  the  spin  exchange  pro- 
cess. The  series  of  cells  contain  Rb®^  metal,  0.05  Torr  Xe-^9 

and  Nj  from  4.3  Torr  to  80  Torr.  Results  of  the  measurements 

2 

are  displayed  in  Figure  2-6  where  we  have  plotted  0 V versus 

GX  rGi 

pressure.  The  bars  on  the  data  points  represent  a 15  percent 

2 

uncertainty  both  in  Nj  pressure  and  °exvre^>  and  are  intended 
only  as  a gauge  of  the  validity  of  the  Least  Squares  fit,  solid 
line,  of  the  data  points  to  Eq.  (28).  Although  the  agreement 
between  our  model  and  experimental  results  in  this  preliminary 
measurement  is  not  as  good  as  would  be  needed  in  order  to  reach 
definitive  conclusions,  we  are  at  this  time  able  to  make  some 
basic  observations: 

a.  The  spin  exchange  process  between  the  alkali  valence 

electron  and  the  noble  gas  nucleus  is  a function  of  the 
buffer  gas  density  in  the  cell.  This  fact  is  con- 
sistent with  the  model  in  which  the  Van  der  Waals 
molecules  have  a significant  contribution  to  the 
exchange  process. 
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b.  We  obtained  three  parameters  related  to  the  three-body 
collision  rate,  the  interaction  strength  times  the 
molecular  lifetime,  and  the  binary  spin  exchange  cross 
section  from  a Least  Squares  fit  of  our  data  to  Eq.  (28): 

(1)  We  obtained  the  three-body  collision  rate  of  Rb, 

Xe  and  N2  per  Xe  atom.  Evaluating  this  rate  at 
60°C  and  1 Torr  n2  we  find: 

T"1  = 4.5  x 10-4  sec-1  (34) 

2 

Major  argued,  based  on  a hard  sphere  collision 
model,  that  an  upper  bound  for  the  termolecular 
collision  rate  per  Xe  atom  at  6 0°C  per  Torr  of 

N2  iS: 

T”1  = 2.5  x 10"3  sec-1  (35) 

Our  estimate  of  the  molecular  formation  rate  seems 
to  be  in  reasonable  agreement  with  that  deduced 
by  Major. 

(2)  The  average  strength  of  the  interaction  multiplied 
by  the  duration  of  the  interaction  per  Torr  N2  we 
found  to  be: 


(y  y t ~ io 


(36) 


4 


l1 
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If  we  assume  that  the  lifetime  of  the  molecular 
state  is  limited  by  Nj  collisions,  we  can  use  the 
same  gas  kinetic  arguments  as  those  employed  in 
(1)  to  estimate  T; 

t = 0.5  x 10~7  sec/'Torr  Nj  (37) 

Substituting  this  value  into  Eq.  (36)  we  obtain  an 
estimate  for  the  average  strength  of  the  inter- 
action during  the  molecular  lifetime: 

<(y)>  ~ 2 x 10®  sec  ^ (38) 

The  strength  is  in  reasonable  agreement  with  what 
had  been  previously  estimated  for  the  binary 
collisions . ’ ® 

(3)  The  data  fit  yields  a spin  exchange  cross  section 

129 

for  the  Rb-Xe  binary  collisions: 

c = 1.5  x 10  ^ cm^  (39) 

exB 

This  value  is  in  excellent  agreement  with  cross 

129 

sections  measured  for  Rb-Xe  in  cells  buffered 
with  500  Torr  He,  where  one  would  expect  binary 
collision  contributions  only.^ 
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We  have  attributed  the  scatter  in  our  experimental  data  to  two 
main  sources  of  difficulty: 

a.  It  was  discovered  in  the  process  of  filling  this  present 
series  of  cells  that  our  buffer  gas  filling  procedure 
could  result  in  very  large  uncertainties  in  the  fill 
pressure  of  the  buffer  gas.  A new  procedure  has  been 
devised  and  a second  series  of  cells  is  being  completed 
for  which  these  present  measurements  will  be  repeated. 

b.  The  non-spin-exchange  part  of  the  relaxation  rate 

I OQ 

for  Xe  is  mainly  due  to  the  interaction  of  the  dipole 
moment  of  the  Xe  nucleus  with  the  magnetic  field  gradi- 
ents present  in  the  cell.  This  relaxation  is  of  the 
form:20 


T 1 = y2D  G2t 2 
2G 


where  D is  the  diffusion  coefficient  of  the  Xe  through 
the  buffer  gas,  G is  the  magnetic  field  gradient  along 

the  quantization  axis  and  t is  the  diffusion  time  from 

wall  to  wall  of  the  atom  in  the  cell.  We  approximate 

- 21 
Tv  by: 

Tv  = R2/ttD  (41) 

where  R is  the  radius  of  the  cell  and  we  take  the  pres- 

2 2 

sure  and  temperature  of  the  diffusion  coefficient  as 


i 

l 
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D P 

D = -2-2  (t/T  ) 
p ' o 


where  p is  the  pressure  of  the  buffer  gas,  pQ  is  stan- 
dard pressure,  T is  the  temperature  of  the  cell,  Tq  is 
273°K  and  a is  a parameter  between  1.75  and  2.  The 


spin  relaxation  due  to  magnetic  field  gradients  can 


then  be  written  as: 


— 1 O' 

T = K P/T 

Cj 


where  K is  some  constant.  From  Eg.  (43)  one  can  see 


that  in  very  highly  buffered  cells,  relaxation  due  to 


magnetic  field  gradients  can  dominate  the  spin  exchange 


process  and  thus  cause  large  errors  in  the  spin  exchange 


measurements , 


To  determine  whether  magnetic  field  gradient  relaxation  was  a 


significant  part  of  the  total  observed  relaxation  in  our  cells, 


we  plotted  the  rate  intercept  found  in  extracting  the  spin 
exchange  information  (see  Figure  2-5)  versus  the  Nj  pressure, 
as  shown  in  Figure  2-7.  The  solid  line  in  Figure  2-7  represents 


a Least  Squares  fit  of  the  data  points  as  a function  of  ^ pres- 


sure according  to  Eq.  (43).  From  the  slope  of  the  line,  we  have 


inferred  that  there  is  a magnetic  field  gradient  in  our  apparatus 


of  approximately  4 [iG/ cm.  This  implies  that  our  present  data 
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N2  PRESSURE  (TORR) 

ERROR  BARS  REPRESENT  A 15  PERCENT  UNCERTAINTY  IN  BOTH  THE 
N9  PRESSURE  AND  THE  DECAY  RATE. 
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reduction  technique  probably  loses  validity  around  ^ pressures 
of  100  Torr  due  to  the  significant  temperature  contributions  as 
shown  in  Eq.  (43).  At  the  present  time  we  are  working  to  reduce 
this  magnetic  field  gradient  in  our  apparatus  by  means  of  gra- 
dient compensating  coils. 

2 . 5 SUMMARY 

We  have  developed  a model  for  the  role  of  the  Rb-Xe  Van  der  Waals 
molecular  complexes  in  the  Rb  electronic  Xe  nuclear  spin  exchange 
interaction.  The  preliminary  measurements  of  the  spin  exchange 
parameters  have  revealed  a significant  contribution  to  the  spin 
exchange  process  from  the  formation  of  the  Van  der  Waals  molecules 
however,  we  are  not  able  to  make  a definitive  comparison  between 
the  model  and  the  experimental  results  at  this  time. 

2.6  CONSIDERATIONS  FOR  FUTURE  WORK 

Our  aim  is  to  reduce  those  sources  of  error  in  our  apparatus 
and  measurement  technique  which  have  been  identified  as  a result 
of  this  series  of  measurements.  Similar  measurements  will  be 
made  extending  the  buffer  gas  regime  out  to  several  hundred  Torr 
in  order  to  check  the  limiting  regions  of  our  model. 

Future  work  will  also  involve  developing  an  understanding  of  the 
role  of  both  the  alkali  nuclear  spin  and  the  noble  gas  spin  by 
considering  similar  measurements  employing  both  Rb®®  and  Xe’*’®'*-. 
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SECTION  III 

NOBLE  GAS  RELAXATION  EFFECTS 

3 . 1 INTRODUCTION 

Recent  measurements  in  our  laboratory  of  the  transverse  decay 
83 

rate  of  Kr  showed  that  the  decay  of  the  nuclear  polarization 

is  dependent  on  the  particular  orientation  of  the  NMR  cell  in 

the  experimental  apparatus.  This  effect  was  previously 
1 2 201 

observed  ’ in  the  relaxation  of  Hg  . The  particular  orien- 
201 

tation  of  the  Hg  cell  which  resulted  in  the  longest  T ^ decay 

time  of  the  Hg  polarization  was  termed  the  "magic  angle"  where 

the  angle  was  between  the  cell  axis  and  the  externally  applied 

magnetic  field.'*'  No  such  angular  effects  have  been  observed  in 
129  199 

either  Xe  or  Hg  both  of  which  have  nuclear  spin  I = 1/2, 

this  suggests  that  the  angular  effect  results  from  quadrupolar 

83  2 01 

type  interactions  present  in  Kr  and  Hg  decay.  In  this 
section  we  will  consider  possible  relaxation  phenomena  associated 
with  quadrupolar  interactions  present  in  nuclei  with  spin  greater 
than  1/2,  developing  particular  expressions  for  nuclear  spin 
species  I = 9/2,  Kr®^ , and  I = 3/2,  Xe'*'^'*'  and  Hg^^'*'. 

3 . 2 THE  EXPERIMENTAL  OBSERVATIONS 

83 

Typical  precessional  decay  data  for  a Kr  sample  is  displayed 
in  Figures  3-1  and  3-2.  This  particular  sample  contained 

QO 

natural  Rb  metal,  150  Torr  Nj  and  5 Torr  Kr  in  a 1 ml  Corning 
1720  glass  cell.  Figures  3-1  and  3-2  represent  data  taker,  at 
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the  same  temperature,  approximately  60°C,  after  an  initial 
constant  pump  time  with  each  plot  corresponding  to  a different 
orientation  of  the  cell  in  the  experimental  apparatus.  The 
data  stream  in  each  plot  is  500  seconds  long.  The  angle  indi- 
cated on  each  plot  is  taken  between  the  cell  axis  and  the 
externally  applied  magnetic  field.  The  cell  axis  was  taken  to 
be  along  a line  that  passes  symmetrically  through  the  tip-off 
region  of  the  cell  as  a matter  of  convenience. 


3.2.1  General  Features  of  the  Data 

We  enumerate  general  observations  concerning  the  data: 

a.  T2  decay  of  the  Kr®^  spin  system  is  a function  of 
angular  orientation  of  the  cell  in  the  experimental 
apparatus. 

b.  The  Kr  signals  cannot  be  solely  described  by  an 
exponential  decay.  In  addition,  the  non-exponential 
character  of  the  decay  is  a function  of  angle. 

c.  The  precession  frequency  is  a function  of  time  and 
cell  orientation.  This  can  be  seen  by  comparing 

{ } 

Figures  3-la  and  3-lb.  Also  comparing  Figures  3-la 
and  3-lc  one  sees  an  indication  of  a possible  phase 
reversal  in  the  second  cycle  of  Figure  3-lc. 

d.  The  decay  characteristics  at  a particular  angle  designa- 
tion are  a function  of  the  position  of  the  Rb  metal 
reservoir  in  the  cell.  In  the  sample  shown,  the  magic 
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angle  for  the  cell  was  at  or  near  92  degrees.  The 
magic  angle  for  the  cell  can  therefore  be  changed  by 
moving  the  Rb  distribution  on  the  cell  walls. 

e.  Previous  measurements  of  the  Kr  decay  characteristics 
as  a function  of  cell  size  demonstrated  a significant 
contribution  to  the  relaxation  processes  from  wall 
effects . 

f.  The  initial  amplitudes  appear  to  be  constant  as  a 
function  of  angle,  implying  that  the  angular  effect 
does  not  significantly  contribute  to  the  T^  processes 
of  the  Kr  spin  system. 


3.3  TRANSVERSE  RELAXATION  EFFECTS 

Transverse  (i.e.,  T2)  relaxation  refers  to  the  loss  of  phase 

coherence  in  a precessing  spin  system.  Spin  de-phasing  in  a 

quadrupolar  interaction  can  arise  through  both  an  inhomogeneous 

broadening  of  the  nuclear  Zeeman  levels  and  transitions  among 

3 

the  energy  levels. 


3.3.1  Level  Shifts  Due  to  a Quadrupolar  Interaction 
We  take  the  full  spin  hamiltonian  to  be: 


H Hm  + Hq 


where  HM  is  the  static  magnetic  interaction: 


hm  = -Yh»0  h 
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with  y the  nuclear  gyromagnetic  ratio,  Hq  the  external  magnetic 
field,  taken  to  be  along  the  z-axis,  and  I the  z-component  of 
nuclear  spin.  Hq,  the  quadrupole  hamiltonian,  represents  the 
interaction  of  the  nuclear  quadrupole  moment  with  electric  field 


gradients.  It  can  be  written: 


4,5 


eQ ^ a2  V 

*Q  ~ 61(21-1)  9x^Sx, 


= 


j.k 


D * 


[3/2(1  . 


\ + Vj1  - 6jkI(I+1). 


(3) 


The  indices  j and  k refer  to  axes  x,  y and  z.  Q is  the  quad- 
rupole moment  of  the  nucleus,  e is  the  elemental  charge, 

2 

S V/ax  .ax,  represents  a component  of  the  electric  field  gradient 
1 

tensor.  To  simplify  the  calculation,  we  limit  ourselves  to 


electric  field  gradients  with  cylindrical  symmetry,  and  desig- 
2 2 

nate  eq  = a V3z'  , where  z'  is  then  the  cylindrical  axis  of  the 
field  gradient.  The  quadrupole  hamiltonian  becomes: 


H„  = 


2 | 

4I(2I?1)  W2  <3  cos2«  - 11  [3Iz  - 1(1  + 3>] 

+ 3/2  sin0  cos 0 [IZ(I+  + I_)  + (I+  + I~)IZJ 

2 2 ) 

+ 3/4  sin20  (I+  + I”  )! 


(4) 


ra 
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where  0 is  the  angle  between  the  z and  z'  axes,  and  I are  the 
raising  and  lowering  operators  associated  with  the  spin  operator 

fa 

I,  in  the  cartesian  reference: 


I-  = I ± i I 
x y 


(5) 


With  the  following  definitions: 


a = (3  cos  0-1) 


b = sin  0 cos  0 


• 2 

c = sin  0 


(6a) 
(6b) 
( 6c) 


a = 


2 _ 
e Qq 


41(21-1) 


and 


= -H 


( 6d) 

( 6e ) 


the  hamiltonian  matrices  for  the  spin  systems  I = 3/2  and 
I = 9/2  are  as  shown  in  Figures  3-3  and  3-4,  respectively. 


The  energy  levels  of  the  system  are  found  through  second  order 
perturbation  theory:7 

E(2)  = E°  + <m  l H I m>  + £ ^ ^H-  (7) 


E - E . 
m 3 
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Figure  3-3.  Hamiltonian  Matrix  for  Spin  I = 3/2  for 
Hamiltonian  Defined  by  Eq.  (1) 

Basis  States  are  Taken  to  be  |l,  mj),  with  Terms  Defined 

in  Eqs.  (6a)  to  (6e) 
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where  we  have  taken  the  quadrupole  interaction  as  a perturbation 
on  the  Zeeman  term.  From  Eq.  (7)  with  the  aid  of  Figures  3-3 
and  3-4  we  can  construct  an  energy  level  diagram,  as  shown  in 
Figures  3-5  and  3-6  for  the  spin  system  1=3/2  and  I = 9/2, 
respectively. 

3.3.2  Effects  of  Energy  Level  Shifts  on  the  Precessinq  Spin 
System 

We  experimentally  observe  the  spin  de-phasing  of  a precessing 
nuclear  sample  by  first  creating  a longitudinal  nuclear  polari- 
zation through  spin  exchange  with  optically  pumped  alkali  vapor 
and  then  rotating  the  external  magnetic  field  by  90  degrees  to 
initiate  precession  of  the  nuclear  spins.  We  equivalently 
describe  this  process  in  terms  of  the  noble  gas  density  matrix. 
We  restrict  ourselves  here  to  nuclear  spin  I = 3/2;  however, 
analogous  equations  can  be  written  for  the  spin  I = 9/2  case. 


We  assume  that  as  a result  of  the  alkali  noble  gas  spin  exchange 
interaction,  the  noble  gas  nucleus  attains  some  low  degree  of 
polarization.  In  the  low  polarization  limit,  the  density  matrix 

Q 

for  the  noble  gas  system  can  be  written: 


P = 


1/4  + 3/2?  0 0 

0 1/4  + 1/2?  0 

0 0 1/4  - 1/2? 

0 0 0 


0 

0 

0 

1/4  - 3/2? 


(8) 
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1ST  ORDER  2ND  ORDER 

ZEEMAN  LEVELS  QUADRUPOLE  CONTRIBUTIONS 


ENERGY  LEVELS  DEFINED  BY  EQ.  (7)  UNDER  THE 
HAMILTONIAN  DEFINED  IN  EQ.  (1)  WITH  TERMS 
DEFINED  IN  EOS,  (6a)  TO  (6e). 


Figure  3-5.  Energy  Level  Diagram  1=3/2 
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ENERGY  LEVELS  DEFINED  BY  EQ.  17)  UNDER  HAMILTONIAN 
DEFINED  IN  EQ.  (1).  WITH  TERMS  DEFINED  IN  EOS.  (6a)  TO  (6e). 


Figure  3-6.  Energy  Level  Diagram  I = 9/2 
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where  J is  a small  number  compared  to  unity.  The  factor,  1/4, 
is  the  equilibrium  population  density  of  each  level.  We  can 
rewrite  the  initially  prepared  density  matrix: 


3/2  0 0 0 1 

0 1/2  0 0 

P = 1/4  [i]  + t ( 9 ) 

0 0 -1/2  0 

00  0 -3/2 


where  [I]  is  the  identity  matrix. 

Rotating  the  magnetic  fields  in  the  experimental  apparatus  is 
equivalent  to  a rotation  of  the  density  matrix,  which  is  given 
by: 

P'  = Dp  Dfc  (10) 

g 

where  D is  a Wigner  D-matrix.  The  D-matrix  for  a rotation  of  a 
spin  3/2  system  about  the  y-axis  through  an  angle  p is  shown  in 
Figure  3-7.  Substituting  Eq.  (9)  into  Eq.  (10)  and  using  Fig- 
ure 3-7  we  find  the  initial  state  of  the  rotated  density  matrix: 

cos  (3  - \/~3  sin  j3  0 0 

-v/T  sin  fi  cos  (3  -2  sin  0 0 1 

P'  (0)  = 1/4  [l]+  f 

-2  cos  (3  -2  cos  /3  - y/T  sin  (3 


0 


Figure  3-7.  Wigner  D-Matrix  for  Rotations  About  the  Y-Axis  of 
Right-Handed  Cartesian  System,  Through  an  Angle p 
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The  precessional  characteristics  of  the  spin  sample  are  described 
by  the  time  evolution  of  the  density  matrix.  The  time  dependence 
of  the  density  matrix  is  given  by: 

p'  (t)  = U(t,  0)  p'  (0)  Ufc(t,  0)  (12) 

where  the  evolution  operator,  U(t,  0)  is  defined  through  the 
differential  equation: 

U(t,  0)  = -i  U(t,  0 ) H (13) 

with  the  condition: 

U(  0 , 0)  = [i]  (14) 

In  solving  Eq.  (13),  we  have  assumed  that  the  quadrupolar  con- 
tribution is  small  so  that  we  may  maintain  the  original  eigen- 
vectors of  the  magnetic  hamiltonian  while  taking  the  eigenvalues 

% 

through  second  order  perturbation.  The  evolution  operator 
becomes : 


(15) 

The  time  dependence  of  the  density  matrix,  using  Eqs.  (11),  (12) 
and  (15),  is  displayed  in  Figure  3-8. 
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Figure  3-8.  Time  Dependence  of  the  Noble  Gas  Density  Matrix 
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The  signal  observed  in  our  experiments  corresponds  to  the 
expectation  value  of  the  x-component  of  the  noble  gas  polariza- 
tion, <IX>>  defined  to  be: 

<V  = Tr[!x  P']  (16) 


with  the  operator  I for  the  spin  3/2  system  written  in  matrix 
form  as: 


1 = T 

x 2 
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(17' 


Thus  our  experiment  is  found  to  correspond  to  the  following 


equation: 

<ix>(t>  = 

_ sin  13  | 

^ 2 1 

cos  (E3^2  “ El/2^  + ^ cos 

(El/2  " E-l/2)fc 

+ 3 cos  ( E_1/2  " E_3/2)t 

(18 

For  the  case  in  which  the  guadrupole  interaction  is  not  present 
the  energy  differences  all  become  the  same  and  Eg.  (18)  col- 
lapses to:  v 


<Ix>(t)  = -5  f sin  13  cos  u>t 


(19 


which  simply  describes  free  precession  of  the  spin  system  in  a 
magnetic  field. 
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Writing  out  the  frequency  terms  in  Eq.  (18)  we  have: 


„ 54  »2b2 

E3/2  " El/2  r“+3aa  + 


(20a) 


2 2 2 2 

E E - Cl’  54  a h I 27  Q ° 

El/2  E-l/2  w 4w 


(20b) 


. 54  u2b2 

E-l/2  " E-3/2  = W ' 3 “a  + 1 


(20c) 


noting  that  the  terms  a,  b and  c have  angular  dependence.  From 
Eq.  (18),  we  see  that  a quadrupole  type  interaction  results  in 
three  precession  frequencies  for  a spin  3/2  system.  In  general 
then,  a spin  system  described  by  spin  I will  have  21  precession 
frequencies  in  the  presence  of  a quadrupole  type  interaction. 

The  different  frequencies  can  be  seen  to  be  not  only  a function 
of  the  strength  of  the  quadrupole  interaction  but  also  of  the 
relative  orientation  of  the  electric  field  gradient  with  respect 
to  the  external  magnetic  field  direction. 

Examples  of  the  precession  signal  described  by  Eq.  (18)  are 
displayed  in  Figures  3-9  and  3-10.  Figure  3-9a  is  an  example 
of  the  free  precession  signal  in  the  absence  of  any  quadrupolar 
interaction,  presented  for  the  purpose  of  comparing  the  effects 
of  the  presence  of  the  quadrupolar  interaction.  The  strength 
of  the  quadrupole  interaction  is  the  same  in  Figures  3-9b  to 


3-18 


Figure  3-9.  Free  Precession  Signal  of  a Spin  3/2  System  in  the  Presence 

of  a Quadrupolar  Interaction 


Precession  Signal  of  a Spin  3/2  System  in  the  Presence 
of  a Quadrupolar  Interaction 
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3-10C  and  was  taken  to  be  2 percent  of  the  magnetic  interaction. 
The  plot  displays  the  effect  of  the  angular  orientation. 

Figure  3-9b  was  attained  by  forcing  the  first  order  contribution 

of  the  quadrupolar  interaction  to  zero.  This  is  accomplished 

2 

by  setting  the  term  (3  cos  9 - 1)  to  zero.  The  angle  d is  found 
to  be  approximately  55  degrees,  which  has  previously  been  desig- 
nated as  the  'magic  angle'.1  Figure  3-9c  corresponds  to  an 
angular  orientation  of  75  degrees,  about  20  degrees  off  the  magic 
angle;  and  Figures  3-10a,  3-10b  and  3-10c  represent  orientations 
of  80,  85  and  90  degrees,  respectively.  The  tick  marks  on  all 
the  plots  designate  the  time  position  of  the  positive  peaks  of 
the  free  maqnetic  precession  signal. 

83 

Various  aspects  of  the  Kr  data  can  now  be  explained  in  terms 
of  our  model,  albeit  the  model  was  derived  for  a spin  3/2 
system; 

a.  There  is  an  apparent  decay  of  the  precession  signal  as 
a function  of  angle.  It  is  apparent  in  the  sense  that 
Eq.  (18)  describes  three  oscillatory  signals  which  beat 
with  one  another.  The  beat  envelope  then  is  inter- 
preted as  a decay.  If  we  were  to  extend  the  time  axis 
on  Figures  3-9c  to  3-10c  one  would  see  the  signal 
increase  again  in  time.  This  phenomenon  is  not  seen 
in  the  actual  Kr  data  due  to  the  presence  of  other 
decay  processes  which  drive  the  signal  amplitude  to  zero. 
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b.  The  beat  envelope  is  not  an  exponential  function.  The 
Kr  data  is  some  combination  of  exponential  decay  due 
to  'real'  decay  processes  modulated  by  the  beat  fre- 
quency. Near  the  magic  angle,  the  signal  is  well 
described  by  an  exponential  function  and  quadrupolar 
effects  are  a minimum,  while  at  angles  more  distant 
from  the  magic  angle  the  decay  signal  develops  a more 
non-exponential  character. 

c.  Using  the  tick  marks  on  the  plots  as  a reference  one 
can  see  that  the  precession  frequency  becomes  shorter 
with  time  as  a function  of  angle.  Figures  3-10a 
through  3-10c  show  a phase  reversal  in  the  decay  signal. 

3.3.3  Origin  of  the  Angular  Effects 

We  noted  earlier  that  the  Kr  decay  data  was  influenced  by  the 
position  of  the  excess  Rb  metal  on  the  walls  of  the  cell  and 
that  the  Kr  spin  relaxation  rate  was  a function  of  the  cell 
size.  We  understand  these  facts  to  imply  that  the  Kr  experi- 
ences a quadrupole  interaction  in  collisions  with  the  cell  walls. 

It  has  been  argued  that  in  a collision  of  a noble  gas  atom  with 

a cell  wall,  the  atom  is  temporarily  bound  to  the  wall  through 

3 

a Van  der  Waals  type  interaction.  The  depth  of  the  well 
between  the  atom  and  a particular  wall  site  is  proportional  to 
the  surface  adhesion  energy.  The  electronic  charge  cloud  of  the 
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atom  undergoes  a distortion  during  the  atom-wall  interaction 
which  results  in  an  electric  field  gradient  at  the  nucleus, 
proportional  to  the  curvature  of  the  potential  well. 

We  can  divide  the  binding  sites  on  the  cell  walls  into  two 
distinct  classes;  one  in  which  the  noble  gas  atoms  bind  directly 
to  the  glass  wall  and  the  other  in  which  the  noble  gas  atoms 
bind  to  the  layer  of  deposited  Rb  metal.  One  would  expect  that 
the  adhesion  energies  between  noble  gas  atoms  and  glass  to  be 
different  from  those  of  noble  gas  atoms  and  alkali  metal.  Since 
it  is  observed  that  the  alkali  metal  does  not  uniformly  distrib- 
ute over  the  glass  surface  then  it  must  be  true  that  there  is 
an  angular  distribution  of  adhesion  sites  on  the  cell  wall. 

Notice  that  this  argument  does  not  depend  on  the  strength  of  the 
binding  energy  of  the  noble  gas  atom  glass  sites  compared  to  the 
noble  gas  atom  alkali  metal  sites;  only  that  the  adhesion  be  dif- 
ferent between  the  two  classes  of  sites.  In  addition,  the  angu- 
lar distribution  of  binding  sites  will  be  only  weakly  dependent 
on  the  macroscopic  shape  of  the  cell.  The  angular  dependence 
observed  in  the  Kr  data  results  from  the  averaging  over  the 
binding  sites  with  respect  to  the  z-axis  defined  by  the  external 
magnetic  field;  as  the  cell  is  rotated,  the  averaging  with  respect 
to  the  space  stationary  z-axis  changes  and  yields  the  observed 
angular  effect.  Redistributing  the  excess  Rb  metal  on  the  cell 
wall  changes  the  angular  distribution  of  the  binding  sites  and 
thus  changes  the  angular  effect  at  a particular  cell  orientation. 
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Wall  collision  effects  on  spin  systems  exhibit  a strong  dependence 
on  the  surface  to  volume  ratio  of  the  cell.  This  general  feature 
has  been  observed  in  our  experiments,  which  gives  us  confidence 
in  our  model  of  the  Kr  wall  collisions. 

3.3.4  Level  Transition  Effects 

The  true  relaxation  rate  of  a spin  system  through  a quadrupole 
coupling  of  the  nucleus  with  the  electric  field  gradients  which 
appear  at  the  nucleus  during  the  interaction  of  an  atom  on  a 
wall  has  been  worked  out  for  Hg^0^  by  Cohen-Tannoud ji . ^ The 
time  evolution  of  the  nuclear  density  matrix  is  given  by  the 
formula: 


§£  = - f * (t ) , [hq(  t-  r)  , p(t)]j  dr  (21) 

•'o 

which  is  second  order  in  the  quadrupole  hamiltonian  as  compared 
to  the  pseudo-relaxation  through  level  shifting  which  is  first 
order  in  the  interaction  hamiltonian. 

Eq.  (21)  also  governs  the  longitudinal  effects  due  to  the  quad- 
rupole interaction.  The  degree  of  angular  effect  resulting  from 
relaxation  described  by  Eq.  (21)  would  be  readily  apparent  in 
the  measured  T^  rate;  however,  no  angular  effects  have  been 


observed  in  the  Kr  Tn  measurements  and  thus  we  conclude  that  we 


(23b) 


LU  GUIDANCE  & CONTROL  SYSTEMS 

Litton  5500  Canoga  Avanue,  Woodland  Hills.  California  01365 


FSCM  06481 


and  thus  the  quadrupoie  relaxation  described  by  Eq.  (21)  is 
characterized  by  a single  rate  constant,  J,  in  our  experiments. 
From  Eq.  (13b)  of  Ref.  5,  the  decay  of  the  transverse  polariza- 
tion can  be  shown  to  be  given  by  the  following  formula: 

<i  > = -2  J <1  > (25: 


which  is  a simple  exponential  decay  function.  The  transverse 
decay  of  a spin  system  due  to  a quadrupoie  interaction  is 
described  by  the  function  represented  in  Eq.  (18)  multiplied 
by  an  exponential  decay. 

3 . 4 SUMMARY 

We  have  experimentally  measured  the  pumping  and  decay  character- 
istics  of  Kr  as  a function  of  the  orientation  of  the  Kr  cell 
in  our  experimental  apparatus.  We  have  qualitatively  explained 
our  observations  based  on  a model  in  which  the  nuclear  spin  is 
3/2.  In  addition,  we  have  written  the  analogous  equations  for 
spin  9/2;  however,  due  to  the  complexity  of  the  spin  9/2  system, 
we  have  not  been  able  to  stabilize  our  computer  codes  enough  to 
give  us  reliable  Least  Squares  fits  to  our  Kr  data  based  on  our 
model.  We  do  not  feel  that  this  is  due  to  any  inherent  diffi- 
culty in  our  model. 
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3 . 5 CONSIDERATIONS  FOR  FUTURE  WORK 

We  plan  to  measure  the  decay  characteristics  of  Xe^-*-  in  a 

83 

manner  similar  to  that  done  for  Kr  . Due  to  the  inherently 

simpler  properties  of  the  spin  3/2  system  as  compared  to  that 

of  the  spin  9/2  system,  we  feel  that  we  will  be  able  to  perform 

131 

parametric  fits  to  the  Xe  data  and  realize  a quantitative 
verification  of  our  relaxation  model. 
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